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Abstract 

CFD has played a major role in the resurgence 
of hypersonic flight, on the premise that numeri- 
cal methods will allow us to perform simulations 
at conditions for which no ground test capability 
exists. Validation of CFD methods is being estab- 
lished using the experimental data base available, 
which is below Mach 8. It is important, however, 
to realize the limitations involved in the extrapo- 
lation process as well as the deficiencies that 
exist in numerical methods at the present time. 
Current features of CFD codes are examined for 
application to propulsion system components. The 
shortcomings in simulation and modeling are identi- 
fied and di scussed . 

Introduction 

An overwhelming degree of reliance has been 
placed on computational fluid dynamics in the 
achievement of hypersonic flight with a NASP-like 
vehicle (National Aero-Space Plane Program). This 
reliance covers the range of design activities; 
from the design of the aircraft configuration to 
the cjesign of the integrated engine system. The 
belief that CFD can be used to predict all of the 
relevant flow physics and chemistry, from aircraft 
takeoff to orbital speeds and return, has been one 
of the principal reasons for the resurgence of 
hypersonic research. 1 Computational methods do, in 
fact, provide us with the unique ability to perform 
ground simulations at high Mach numbers for which 
no ground test capability exists. Above Mach num- 
bers of approximately 8, ground test facilities do 
not duplicate the relevant flight simulation param- 
eter such as Mach number, Reynolds numbers, gas 
composition, and enthalpy level. Numerical analy- 
sis remains as the principal approach to the design 
of the aircraft and the propulsion system. It is 
possible that some data could be obtained from 
rocket test vehicles. These data would be limited, 
however to such items as the state of the boundary 
layer, boundary layer transition location, length 
of the transition zone and surface heat transfer. 
Testing of a scramjet propulsion system on a rocket 
vehicle would present a major problem due to the 
fact that scaling of the combustion process is not 
feasible. Therefore, a full sized propulsion mod- 
ule would be required as the test article on a 
rocket test vehicle. Such testing might be better 
approached through use of larger vehicles such 
as the Space Shuttle. However, it is currently 
believed that the costs associated with the flight 
testing described above would be extremely high. 

CFD, therefore, remains as a viable alternative. 

It must be pointed out that the philosophy regard- 
ing CFD is based on the fact that an experimental 
data base exists below Mach 8. Those data provide 
the means for assessing the accuracy of the numeri- 
cal methods, as well as for calibrating the codes. 
Extrapolation to flight conditions, where no data 


‘Deputy Chief, Internal Fluid Mechanics 
Division; Associate Fellow, AIAA. 


exists, can then be carried out. In theory, there- 
fore, the procedure described would provide us with 
ground simulation throughout the Mach number range 
from takeoff to orbital velocities. It is crucial, 
however, to realize the capabilities and the limi- 
tations involved in the extrapolation process as 
well as the numerical methods at the present time. 
In this paper, the current features of the numeri- 
cal methods available for analyzing the flow in 
multi-element hypersonic propulsion systems will 
be presented. Since this paper has as its focus 
the propulsion system, the aircraft configuration 
is only discussed relative to the propulsion flow 
path. Only the integration effects of the fore- 
body and nozzle afterbody configuration will be 
addressed. Emphasis is placed on the inlet behav- 
ior, combustor and nozzle characteristics includ- 
ing ionization and dissociation effects, as well 
as finite rate and equilibrium chemistry effects. 
The critical limitations introduced by the need for 
turbulence, boundary layer transition and chemical 
modeling are discussed; and the effect of various 
models is illustrated for the propulsion system 
components. Shortcomings associated with the 
extrapolation to higher speeds are also presented. 
Current and future activities, which are directed 
toward improving the modeling, are discussed. 

These activities include elements such as direct 
numerical simulation, wall and shear layer turbu- 
lence modeling and probability density functions 
for reacting flows. 

Propulsion CFD Validation 
Hypersonic Propulsion System 

In order to demonstrate the ability of CFD 
codes to perform propulsion system computations and 
to examine their strengths and their shortcomings, 
a typical propulsion configuration must be chosen. 
In this paper, a combined ramjet/scramjet system 
is assumed which has a common flow path. The low 
speed or ‘'accelerator" portion of the engine, which 
would provide sufficient speed for ramjet opera- 
tion, is not considered here. The assumed flight 
conditions would correspond to subsonic combustion 
ramjet operation from flight numbers of approxi- 
mately Mach 3 up to approximately Mach 5.5, where 
supersonic combustion operation would commence. 
Supersonic flow within the engine results in lower 
temperatures and pressures, thereby reducing heat 
fluxes and internal forces. The lower temperatures 
also allow heat to be added from the fuel without 
imposing high air dissociation losses. However, 
ignition and combustion of the fuel within the com- 
bustor remains a major challenge. The idea of an 
oblique detonation wave ramjet has been proposed, 
wherein fuel is injected upstream (inlet) for pre- 
mixing with air. The mixture is then ignited by a 
shock wave within the combustor. Detonation wave 
stability, completeness of combustion and premix- 
ing feasibility issues must be solved to under- 
stand this concept. In this paper, the detonation 
engine computations are not considered. Our inves- 
tigation, therefore, centers on the use of CFD for 
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an integrated airbreathing propulsion system oper- 
ating at altitudes from near earth to the limits 
of the stratosphere. In the absence of specific 
ramjet/scramjet designs, a generic type propulsion 
system will be used which includes many of the 
features which are considered important in real 
systems. Those engine features which cannot be 
analyzed, for example, complex combustor geome- 
tries, will be discussed as limitations of CFD for 
hypersonic propulsion systems. 

The level of sophistication required in the 
numerical methods varies for each of the system 
components. For the blunt forebody, a thin layer 
Navi er-Stokes code is required at the nose to 
handle the strong viscous shock, followed by para- 
bolized Navier-Stokes method to the inlet. Around 
Mach 7, the vehicle generates a high temperature, 
high enthalpy flow field. The elevated tempera- 
tures in the boundary layers of the vehicle disso- 
ciates the oxygen and nitrogen molecules into their 
respective atomic species which may also become 
ionized. The process can be catalyzed and enhanced 
further by vehicle surface contaminants. This dis- 
sociation process may lead to variations in local 
flow-field molecular weight, specific heat and 
transport properties with resultant associated 
changes in dynamic airflow, heat transfer and 
combustor kinetics. These effects cannot be accu- 
rately modeled in ground test facilities. Prop- 
erties of shock heated air and combustion gases 
are needed for performance predictions. Three- 
dimensional thin layer or Navier-Stokes are 
required for the flow in the inlet with the pres- 
ence of multiple shocks and possible flow separa- 
tion and unsteadiness. The combustor requires a 
Navier-Stokes solution including finite rate chem- 
istry as does the nozzle flow field. It is obvi- 
ous, therefore, that real gas effects, ionization, 
recombination, nonequilibrium effects, and wall 
catalyticity must be considered at various loca- 
tions through the system. 

CFD Validation 

As mentioned in the Introduction, CFD extrapo- 
lation is required for the design of the aircraft/ 
propulsion system at the higher Mach number range 
(M > 8). A critical activity, therefore, is asso- 
ciated with validation/calibration of the numerical 
techniques. To what extent or level of sophistica- 
tion is validation required? Must the numerical 
codes duplicate all of the physics and chemistry 
in the flow? Before answering these questions, 
one must be aware of the manner in which codes 
will be used in the design process. In the case 
of the inlet, the designer is interested in cer- 
tain performance parameters such as the amount of 
mass captured, adiabatic kinetic energy efficiency, 
pressure recovery, heat load, spillage drag and 
exit profile. Ideally, then a designer can com- 
pute, for a variety of possible inlet geometries 
over a range of Mach numbers, all of the parameters 
cited above and, eventually, arrive at a fairly 
efficient inlet design. At this point, however, it 
is necessary to ask to what extent is the designer 
interested in the physical and chemical phenomena 
occurring in the inlet, such as shock/boundary 
layer interactions, secondary/corner flows, mass 
injection, transition, bleed, equilibrium chemis- 
try, flow separation and unsteady flows due to 
pressure oscillations and shock interactions? The 
answer is perhaps quite straightforward. The 
designer is interested in all of the physics and 


chemistry insofar as it affects the performance 
parameters cited earlier. If the physical /cnemi cal 
phenomena do not affect the inlet performance or 
loadings, it is safe to say that the designer wou ] d 
not be particularly upset if such features were 
missing in the numerical code. 

The relevance of the preceding discussion to 
the issue of code validation/calibration is that 
the extent of the validation process and the number 
of effects to be accounted for are strongly depend- 
ent on the magnitude of the individual effects on 
inlet performance. Another way of stating it is 
to ask which of the items in the list of physical/ 
chemical phenomena must be modeled by the CFD'er 
in order for the designer to have valid answers. 

If none of the items can be ignored, then all of 
them must be included as the numerical codes are 
exercised to develop sensitivities to the indivi- 
dual modeling. One must realize that in addition 
to the modeling required, there are also issues 
related to the numerics and the math modeling. 

These issues involve computational grid sensitiv- 
ity, the ability to capture discontinuities, sensi- 
tivity to internal code parameters, the effect of 
numerical boundary conditions, and conservation of 
mass, momentum, energy and species. Therefore, 
one must be concerned with numer i cal /mathemat i cal 
modeling as well as physi cal /chemi cal modeling. 

If these activities can be properly executed, then 
one may proceed to the next stage, which includes 
the identification of critical experiments and com- 
parison of the numerical results with experimental 
results. Although the experimental methods will 
not be discussed, it is equally important to estab- 
lish the validity and accuracy of the measured 
data. 

The discussion regarding code validation has 
been somewhat general to this point. To be more 
specific. Table 1 shows the critical forebody 
design or performance parameters, the quantities 
that must be computed by the codes and the physics/ 
chemistry modeling requi rements . The same informa- 
tion is shown in Table 2 for the inlet. Table 3 for 
the combustor and Table 4 for the nozzle. From the 
viewpoint of a researcher, it is believed that a 
fundamental understanding of the physics and chem- 
istry within any system must first be understood 
and then modeled. Its relative importance (and 
perhaps some eventual control) must first be under- 
stood, before a massive sensitivity study (grid, 
internal parameters) is carried out for every phys- 
ical or chemical phenomena known to exist in the 
component. It is proposed, therefore, for CFD val- 
idation, that both experimental and numerical 
research should proceed from the basis of develop- 
ing understanding first, secondly, making judgments 
on the importance of various phenomena, and then 
performing numerical sensitivity studies. 

The propulsion system components will be dis- 
cussed sequentially in the following sections of 
the paper. The items of special concern are 
unsteady flow behavior in the inlet, combustion- 
turbulence modeling in the combustor and shear 
layer/boundary layer characteristics in nozzle flow. 

Numerical Methods for Inlet Flow 
Numerical Schemes 

Typical results for a high speed inlet will be 
presented in this section. A fairly significant 
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number of analyses have been carried out over the 
last several years at various laboratories (i.e. 
NASA Lewis, Langley, Rose Engr, APL/JHU). As a 
means of illustrating the current capability of 
inlet codes, we consider the results obtained by 
Benson and co-workers^" 4 using both parabolized 
and Reynolds averaged Navi er-Stokes solvers. The 
PNS code solves the system of equations for hyper- 
sonic or supersonic flow by the linearized block 
implicit scheme of Briley and McDonald. 5 Sijice 
parabolized solvers have inherent limitations 
regarding separation regions, an approximation is 
used to allow the code to march through small 
areas of separated flow. The Reynolds averaged 
Navi er-Stokes code (PARC3D) solves the basic equa- 
tions in strong conservation form with the Beam- 
Warming approximate fac tor i zation algorithm. 6 it 
uses central differencing and Jameson type artifi- 
cial dissipation. 7 Originally developed by Pulliam 
and Steger, 8 ARC was modified for propulsion analy- 
sis by Cooper 9 and accordingly named PARC. The 
simple rectangular inlet configuration shown in 
Fig. 1 was analyzed with the expectation that 
experimental data would be available. A flat 
plate of 30 in. length preceded the entrance to 
the inlet in order to simulate the boundary layer 
growth on the forebody of a hypersonic aircraft. 
Compression wedges form the top and bottom walls 
of the inlet and the contraction ratio was equal 
to 5. Swept sidewalls which connect the upper and 
lower walls prevent compressed flow from spilling 
over the inlet sides. Computations were made at 
an entrance Mach number of 12.25. Various turbu- 
lence models were used in the solutions, including 
those developed by McOonald-Camarata , Bushnell and 
Beckwith and Bal dwi n-Lomax . Two-dimensional PNS 
solutions were carried out using grids of 100 by 
1000, and the three-dimensional cases used 80 by 
60 by 750. The two-dimensional Reynolds averaged 
Navier-Stokes solutions used 100 by 200 grids 
whereas the three-dimensional solutions were car- 
ried out on grids of 150 by 81 by 41. Computations 
were performed on the Cray X-MP and the Cray 2. 

The PNS code described above has been modified by 
Liou 10 to include real gas effects. In addition, 

Yu et al. 11 have incorporated finite rate and 
local equilibrium approaches in the chemical reac- 
tion model for dissociation and ionization of the 
inlet air. The finite rate approach involves the 
simultaneous solution of eleven species equations 
coupled with the fluid dynamics equations. In the 
local equilibrium approach, a chemical equilibrium 
package has been developed and incorporated into 
the flow code to obtain chemical compositions 
directly. Gas properties for the reaction product 
species are calculated by methods of statistical 
mechanics and fit to a polynomial form for spe- 
cific heat . 

Rose and Perkins 1 ^ have used Kumar's explicit, 
time-accurate implementation of MacCormack's algo- 
rithm for solving the full Navi er-Stokes equations. 
Grid sizes used for the computations were 201 in 
the streamwise direction, 61 in the compression 
direction and 27 in the cross stream direction. 

Inlet surfaces were assumed to be nearly adiabatic 
with surface temperatures equal to stagnation 
temperatures . 


Computed Inlet Results 


Contour plots of constant Mach number within 
the inlet, obtained with the PNS analysis are 
shown in Fig. 2. The concentration of lines near 


the walls indicates the boundary layer, while the 
contour concentration in the freestream indicates 
shock wave locations. The flow (right to left) 
features seen in the figure are boundary layer 
buildup on the flat plate followed by thickening on 
the side walls and ramp surface. Shocks generated 
by the compression wedges are seen as horizontal 
lines, and the sidewall shocks as vertical lines. 
The shocks generated by the compression surface 
glance across the sidewall boundary layers, produc- 
ing a thickening of the boundary in the vicinity 
of the shock, and thinning of the layer in the cor- 
ners. Strong secondary flows are developed in the 
shock wave/boundary layer interaction. These 
flows, in turn, affect the downstream inlet flow 
field. A more detailed view of the flow field Mach 
number contours is seen in Fig. 3. In this view, 
which looks upstream toward the inlet entrance, the 
ramp and cowl compression shock waves emanating 
from the leading edges are clearly discernible, as 
well as the secondary flow developed in the cor- 
ners. The boundary layer growth on the wedges and 
sidewalls has been highly distorted by interactions 
with the compression shocks. Benson et al. 4 des- 
cribe their computations in the following way, 
which relates to the important physical processes 
occurring: "The low energy flow of the sidewall 

boundary layer has been swept up the sidewall by 
the ramp shock and down the sidewall by the cowl 
shock. Near the sidewall where the secondary flows 
collide, one sees a secondary shock as a vertical 
line. As the flow proceeds downstream to the cen- 
ter plane, the shock waves cross and are distorted 
by interaction with the sidewall boundary layer and 
the expansion fan on the ramp surface. At this 
station, the ramp shock is reflecting from the cowl 
surface and the cowl shock is seen as the upper 
horizontal line. The lower white horizontal line 
corresponds to the edge of the ramp boundary layer. 
The vortices generated by the shock/boundary layer 
interactions have pulled away from the sidewall 
while interacting with each other. Proceeding to 
the last plane on the left, the expansion generated 
on the lower surface causes a strong pressure gra- 
dient from top to bottom. Low energy flow along 
the sidewall moves into the corner formed by the 
sidewall and the ramp surface. As the shock wave 
created by the ramp and reflected from the cowl 
strikes the ramp surface, the low energy flow in 
the corner separates. The PNS analysis cannot be 
made to proceed farther due to the magnitude of the 
separation in the corner. At this last station, 
the flow is seen to be highly distorted with sepa- 
ration in the lower corners, vortical flow near the 
sidewalls and thick boundary layers on both the 
ramp and cowl surfaces." 

An alternate visualization of the three- 
dimensional flow is obtained with particle tracing, 
as shown in Fig. 4. Computational particles were 
placed along the sidewall of the inlet. Interac- 
tion of the ramp and cowl shocks with the sidewall 
boundary layer causes the particles to converge 
near the shock interaction point. The particles 
are then displaced due to the vortex motion des- 
cribed earlier. Flow migration details are clearly 
depicted in this computational simulation. Since 
this vortex phenomenon persists downstream, it is 
proposed that judicious fuel injector design and 
placement may lead to enhanced mixing of air and 
fuel in the combustor. 

Computations performed by Rose are shown in 
Fig. 5 for a rectangular inlet at Mach 5. The 
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same features shown in the previous inlet are also 
seen in the Mach 5 inlet, i.e., strong viscous 
shock effects leading to large regions of vortex 
flow. Rose carried out a number of numerical 
experiments to control the vortex phenomena in the 
corner regions. Figure 6 shows the baseline or no 
control case, followed by cowl cutback, cowl bleed 
and removal of a part of the sidewall. These modi- 
fications were made near the inlet ramp shoulder. 

It may be seen that these modifications were 
ineffective in eliminating the vortex region. Even 
with the cutback sidewall, the low momentum fluid 
exits along the entire sidewall. Some attenuation 
is seen along the cowl surface for that case. It 
is evident, however, that the shock-boundary layer 
physics within a rectangular shaped inlet will lead 
to pressure losses in the corner regions. However, 
if these regions can be utilized in an “integrated 
design approach," then combustor/nozzle design may 
benefit substantial ly. 

Improvements Required 

At the present time, the design of a low drag 
highly efficient inlet has not been demonstrated. 

It must be noted here that the geometries computed 
are extremely simple and, to a large extent, the 
computations are concentrated on developing a fun- 
damental understanding of glancing shock/boundary 
layer behavior within an inlet. The results 
obtained yield a clear picture of the physics asso- 
ciated with the compression process. Observations 
on high speed inlets verify that the secondary 
flow regions do in fact occur, as predicted by the 
numerical code. So far, therefore, the numerical 
results provide realistic answers. However, real 
inlets are more complex in shape, may have movable 
geometry, will spill flow over a wide Mach number 
range, will probably include bleed at low Mach 
number flight and blowing at the higher speeds. 

In addition, some means for shock control and elim- 
ination of buzz and unstart will be required. A 
design CFD code, then, would require a full three- 
dimensional, time accurate Navi er-Stokes code capa- 
ble of handling spillage and mass addition or 
removal. Currently, these individual effects, such 
as spillage, have been computed for steady inlet 
flows with simplified geometries. 13 However, 
incorporation of all of these features into a time 
accurate computer code will make efficient computa- 
tions all but impossible. It is worth noting that 
the simple geometry shown in Fig. 1 was tested and 
found to exhibit unsteady flow, again stressing 
the need for a time accurate solver. It is fur- 
ther noted that the issue of dispersive and dissi- 
pative errors associated with various time accurate 
finite-difference numerical schemes with high grid 
density is only now being rigorously addressed. 14 
Proper boundary condition treatment is also neces- 
sary to avoid degrading the accuracy of numerical 
solutions. Pure acoustic radiation boundary condi- 
tions do not properly account for the effects of 
unsteady vortex shedding, disturbance amplifica- 
tion caused by separation and the response due to 
both vorticity oscillations and entropy waves. 14 

Up to this point, the discussion has centered 
on the geometric configuration of the inlet. 

There are also a number of problems associated 
with understanding the fundamental flow physics 
and our ability to model them. These physics are 
associated with the transition of the boundary 
layer and the nature of the turbulence modeling. 
Currently, there is no way to determine where the 


boundary layer begins to transition from laminar 
to turbulent, as well as the length of the transi- 
tion zone. The state of the boundary layer at 
hypersonic speeds is unknown, and relami nar i zation 
may occur in an unpredictable manner. Direct 
numerical simulation may provide an approach for 
transition modeling at high speeds. 

Reacting Flow Code 
Numerical Schemes 

A representati ve code for the analysis of com- 
bustor flow fields is described in Ref. 15. The 
three-dimensional code employs an implicit finite 
volume, lower-upper (LU) time marching method to 
solve the complete Navi er-Stokes and species equa- 
tions in a fully-coupled and efficient manner. The 
hydrogen-air chemistry model includes nine species 
and eighteen reaction steps. The code has been 
demonstrated for normal hydrogen injection into 
supersonic airstreams for nonreacting flow and also 
compared to other numerical schemes. The LU code, 
in its two-dimensional version, incorporates com- 
prehensive real gas property models to account for 
high temperature flows, and includes finite rate 
or equilibrium chemistry. The code, RPLUS, is 
formulated based on eigenvalue upwinding. The 
scheme has the efficiency and robustness of an 
implicit scheme, with an operational count compara- 
ble to that of an explicit scheme. This feature 
of the code is of critical importance for three- 
dimensional calculations of a large system of equa- 
tions for reacting flows. Vectori zation of the 
code is performed by a reorganization of the 
indices of the grid points for parallel processing 
planes . 

Now, we consider the case of a simple or crude 
combustor which utilizes some hydrogen injection 
from a single or double wall position into a Mach 4 
freestream (see Fig. 7). For the purpose of per- 
forming computations, we assumed the temperature 
of the hydrogen was set at 700 K and the airstream 
at 1300 K. The boundary conditions include no-slip 
and adiabatic top wall, with gradients of variables 
in the streamwise direction assumed to be zero. 
Symmetry boundary conditions are assumed at both 
side walls and bottom wall. Grid sizes required 
for resolution of the jet interaction features 
were on the order of 60 by 40 by 40 (x,y,z) with 
clustering in both x and z for jet resolution and 
y clustering for boundary layer resolution. It 
should be re-emphasized that the case under consid- 
eration represents only a small portion of a real- 
istic combustor geometry. In other words, this 
rather large computation is only a "unit problem 1 ' 
associated with one of the propulsion components. 

The results obtained with this numerical code and 
the grid discussed yield a great deal of the phys- 
ics inherent in the injection process. Observa- 
tions regarding jet interaction phenomena made over 
the last 20 years are accurately replicated by the 
code. It is noted that in the section entitled CFD 
Validation the ability to capture discontinuities 
was mentioned as one of the important numerics 
issues. In that regard, the recent work by Shuen 
and Liou 1 ^ was directed toward a flux splitting 
algorithm for viscous flows with nonequilibrium 
chemistry. Upwind TVD differencing was used with 
a Roe flux splitting scheme. Nonequilibrium, 
frozen chemistry and idea gas assumptions were used 
in a sample calculation for the single jet injec- 
tion case described above. 
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Computed Combustor Results 

Figures 8 and 9 show a typical Mach number 
and temperature contour on yz planes for various x 
locations. Yu et al. 18 describe their computa- 
tions in the following manner, which relates to 
the important physical /chemical phenomena: "Just 

behind the injecting orifice, the Mach number con- 
tours show a strong bow shock very close to the 
wall. Under the bow shock, the circular Mach num- 
ber contour indicates the existence of the barrel 
shock structures. The jet has been bent and flows 
almost parallel to the primary flow. The penetra- 
tion of the jet increases as the flow goes down- 
stream. The shape of the jet also deforms as the 
flow goes downstream due to the presence of the 
streamwise vorticity in the lee of the injector. 

The secondary motion formed by two counter-rotating 
vortices gives the bent-over jet a kidney shape. 

In the temperature contour plots, the hottest 
region is along the wall because of the viscous 
dissipation in the high speed flow enhanced by the 
combustion of the H 2 and air. Away from the wall, 
by comparing the Mach number plots and the tempera- 
ture plots, it is evident that the temperature 
increases after the bow shock. Further downstream 
a thin flame zone characterized by higher tempera- 
ture is developed. Vigorous chemical reaction 
occurs in the flame zone." Mach number contours 
are also shown for the xy plane across the orifice 
for both the single (Fig. 10) and double jet 
(Fig. 11) injection computation. Features near 
the injection orifice(s) indicated the presence of 
a barrel shock structure. Separation and reattach- 
ment features upstream and downstream of the jets 
are also clearly discernible. In the dual injec- 
tion case, blockage caused by the first jet allows 
the second jet to penetrate further into the free- 
stream due to stronger expansion. This behavior 
is consistent with the penetration correlation pub- 
lished by Povinelli et al., 17 in that the free- 
stream momentum deficit approaching the jet is one 
of the parameters governing jet penetration. 

Improvements Required 

So far, we have shown that for flush, single 
and dual wall injection, numerical computations 
provide realistic answers. The results, however, 
cannot be extended for lengthy downstream dis- 
tances with the grid resolution obtained around 
the jets, due to computer limitations. Hopefully, 
once the jet regions are computed, one could pro- 
ceed with different, coarser grids downstream. 

Some of the general flow field property distribu- 
tions would need to be passed on to the new grid. 

We must further consider that wall jet penetration 
is limited in terms of providing fuel distribution 
in supersonic streams. Flush wall jet penetration 
is on the order of 10 orifice diameters. 17 For 
that reason, devices like struts which span the 
combustor completely 18 0 r partially 19 are required 
for distribution of fuel. The numerical codes 
must, then, be capable of multiple wall and strut 
injection from many points on the solid surfaces. 
With partial swept struts, a typical arrangement 
might include 20 injection ports from the ridqe 
line and trailing edge of each of six struts,' 9 as 
well as selected combustor wall locations of a 
comparable number. Clearly, our numerical ability 
to perform this task is not possible. Other 
approaches to combustor design may be possible. 

For example, the possibility of enhanced mixing by 


the streamwise generation of vorticity has been 
suggested 20 and addressed in a number of ways. 21 * 22 
Vorticity generation has been found successful for 
jet engine exhaust flows 22 and is currently under 
study at many laboratories for supersonic streams. 
Since most of the generation occurs through an 
inviscid mechanism, 23 the concept should work 
equally well in supersonic streams. Since it has 
been shown in the preceding section, Inlets, that 
a typical high speed inlet generates vortex flow, 
this flow phenomenon could be utilized as it 
passes downstream, to enhance combustor mixing. 
Vortex-shock interactions 24 and shock-shear layer 
interactions 28 have also been postulated for the 
purpose of enhanced mixing. However, these phenom- 
ena have so far proven of limited value for the 
improvement of mixing. 

The various phenomena described above are all 
under continued investigation. Research on "explo- 
sive growth" in shear layer stability underway at 
NASA Lewis 28 * 27 may provide some assistance in this 
problem. Regardless of the understanding which 
will develop over the next few years, the numeri- 
cal analysis will require extensive computer capa- 
bilities to handle the geometries which evolve. 

It is unlikely that simple combustor geometries 
will satisfy the requirement for this propulsion 
system. 

So far, we have discussed the combustor geome- 
try to a great extent. There is, however, even a 
greater problem associated with the fundamental 
phenomena of turbulent combustion and our ability 
to model it. Although a great deal of research 
has been carried out, the nature of the interac- 
tion between the fluctuating flow field and the 
chemical reaction steps or processes is not under- 
stood. Simple approaches related to the gross fea- 
tures of relatively ideal burner behavior have been 
made, such as flame propagation in turbulent mix- 
tures. Those features studied have related to 
changes in flame speed due to increasing surface 
area and increasing transport but have not specifi- 
cally dealt with how the flow field modifies the 
chemical reaction scheme. The proper modeling of 
turbulent flow reaction remains as a key require- 
ment for the CFD of combustors. 

For turbulent chemical reacting flows, the 
evaluation of the mean source (or sink) of the 
chemical species due to chemical reactions repre- 
sents a major unsolved difficulty. The formation 
(or destruction) rates are nonlinear functions of 
temperature and species concentrations and thus 
knowledge of the mean-valued properties is insuffi- 
cient to evaluate the mean formation rate. For 
example, in finite rate chemistry models, the mean 
formation rate calculated by using the mean-value 
temperature and species concentration in the Arrhe- 
nius form will lead to errors up to three orders 
of magnitude. Formally, the mean reaction rate 
could be obtained by decomposing the temperature, 
density, and mass fractions appearing therein into 
mean and fluctuating components, then taking the 
time (or Favre) average. The resulting equations 
involve many second order moments which need to 
be solved by extra transport equations or to be 
modeled. In this approach, the neglect of the cor- 
relation greater than second order terms is unsat- 
isfactory and leads to erroneous solutions, while 
the retention of the higher order term renders the 
approach intractable. 28 
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The most convenient way to circumvent this 
problem is to use probability density function 
(pdf) method. About the simplest and most popular 
approach is to specify a two-parameter form of the 
pdf in terms of mean and variance of the conserved 
scalar. The transport equations are readily 
obtained from the conservation of mass and chemi- 
cal species. By solving these two equations, the 
pdf at each location for the whole flowfield can 
be obtained. Next, if the correlation is known 
between the conserved scalar and the chemical spe- 
cies concentra t ion (which is usually calculated by 
equilibrium method), then the distribution of the 
chemical species concentration can be obtained by 
performing a straight forward integration. This 
method strongly couples the turbulence and chemical 
reactions. However, the choice of the pdf form and 
the usage of equilibrium assumption inevitably lead 
to errors in predicted results. There are more 
complex approaches along the same line, for example, 
the two-variable formalism adopted by Janicka and 
Kollmann 29 for the H 2 /a i r flame. A single combined 
reaction progress variable was introduced to des- 
cribe completeness of the three-body reactions while 
the two-body reactions are assumed to be in equi- 
librium. Thus the joint pdf for mixture fraction 
and reaction progress variable is needed to deter- 
mine the thermochemical state of the flowfield. 23 

Nozzle CFD Anal ysi s 
Numerical Schemes 

In order to illustrate nozzle CFD capability, 
a number of numerical schemes may be chosen. 

Ruffin et al. 30 have used a Reynolds averaged 
Navi er-Stokes solver which contains a LDU-ADI 
scheme with Roe averaging and MUSCL differencing. 

The algorithm is diagonal in structure and requires 
minimal CPU per iteration. Laminar nonreacting 
computations of a nozzle exhaust flow field were 
made using two patched grids to model the geometry. 
The first grid consisted of 20 by 99 by 35 grid 
points and the second grid contained 51 by 99 by 
52 grid points. Baysal et al. 31 used two Navier- 
Stokes numerical schemes to compute nozzle flows. 
Both schemes were two-dimensional. One scheme was 
an implicit, upwind solution and constant the 
second scheme was an explicit MacCormack and had 
variable y. Bal dwi n-Lomax turbulence modeling 
was used, and the grid dimensions were 155 by 131. 

Lai and Nelson 32 used the three-dimensional PARC 
Navi er-Stokes code to compute nonaxi symmetri c noz- 
zle flows. The numerical scheme in PARC employs 
three point central differencing uniformly through- 
out the flow field to approximate spatial deriva- 
tives. Second and fourth order Jameson type 
artificial dissipation is included. Diagonaliza- 
tion of the inviscid terms simplified the block 
pentadi agonal system of equations to a scalar pen- 
tadiagonal system. The numerical scheme employs 
an ADI Beam and Warming approximate factorization. 

In general, the boundary conditions assume mid- 
plane spanwise symmetry, with no slip velocity and 
adiabatic wall temperature. For the external far 
field, quiescent air at normal conditions is 
assumed and nozzle exit plane conditions are speci- 
fied. Streamwise flux gradients at the outflow 
boundary are assumed negligible. 

Nozzle Resu 1 ts 

Three-dimensional laminar computations by 
Ruffin et al. 3 ^ were made for a flight Mach number 


of 7.3. The geometry chosen was that of a corre- 
sponding experimental configuration (Fig. 12). 

The grid zones chosen for the numerical study are 
shown in Fig. 13. The first grid is used for cal- 
culations up to the nozzle exit plane followed by 
interpolation of data to the second grid. The 
internal nozzle flow at the exit plane is modeled 
by specifying boundary layer profiles. The compu- 
ted Mach number contours in the symmetry plane and 
the two cross-flow planes are shown in Fig. 14, 
with the forebody at 0°, afterbody ramp in the noz- 
zle block region at 15°, and the long expansion 
plate at 20°. A complex interaction may be seen 
in the plume where the bow shock and the side edge 
vortex come together. Also, the ramp shock gener- 
ated on the windward ramp causes high pressure flow 
to move to the leeward region. Figure 15 shows the 
predicted particle traces on the ramp expansion 
surface. In the words of Ruffin, the results are 
discussed in terms of the physical processes taking 
place: "Figure 15 (original in color) allows the 

identification of several flow features near the 
body surface. The turning of the oil traces corre- 
sponds to turning of the high pressure flow from 
the side edge of the ramp and above the nozzle 
block toward the low pressure expanded flow above 
the ramp. Other traces converge onto the separa- 
tion line of a vortex near the side edge of the 
ramp. Additional traces correspond to a separation 
bubble at the location just after the ramp angle 
transitions from 15° to 20°. Finally, the foot- 
print of the shear layer is shown by the particle 
traces near the far edge of the ramp. It should 
be noted that the predicted separation that occurs 
on the ramp may be induced by the assumption of 
laminar flow. This feature may not be present in 
the experimental flow field, which will be turbu- 
lent." It is noted here that Ruffin's computations 
were carried out in order to provide information 
relative to the design of a nozzle validation 
experiment. The proposed experimental model is 
fairly simple relative to an exhaust module of a 
NASP type vehicle. In spite of the simplicity of 
the model, the computed results very clearly point 
out the complexity of flow patterns that arise 
under these conditions. 

In an effort to provide better definition of 
the nozzle flow field features, use of adaptive 
gridding has been studied by Hsu. 33 Figure 16 
shows the Mach number contours and grid using a 
regular grid (Fig. 16(a) and (b)> and using an 
adaptive grid (Fig. 16(c) and (d)). These computa- 
tions were carried out on a simple nozzle shape 
using a 81 by 201 regular grid. Grid adaption was 
carried out in the y direction only in the region 
above the cowl. Figure 16 shows that the use of 
adaptive grid yields a sharper shock and thinner 
boundary layer as well as diminishing the large 
region of shock induced boundary layer separation 
where the shock intersects the upper wall. 

Baysal 31 also used adaptive gridding to com- 
pute the straight wall nozzle shown in Fig. 17. 
Again, the model is a fairly simple representation 
of scramjet nozzle and afterbody. Results with 
the implicit code are shown in Fig. 18, for nonre- 
acting air from the nozzle and an external Mach 
number of 6. The nozzle/external pressure ratio 
was approximately 5. The shear layer can be seen 
originating at the cowl tip with an expansion, and 
it is deflected upwards at about 13°. Flow expan- 
sion proceeds down the ramp without separation. 
Computed values of the surface pressure on ramp 
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'show good agreement with measured data, Fig. 19. 

The next calculation by Baysal^ 1 was carried out 
with the explicit code mentioned in the section 
Numerical Schemes, which included variable y. The 
exhaust gas used was 50 percent F-12 and 50 percent 
argon; the free stream Mach number was 6. Density 
comparisons are shown in Fig. 20 and pressure con- 
tours in Fig. 19. The expansion process proceeds 
downstream and the shear layer deflects upwards at 
about 15°. The density is higher in the shear 
core, expanding back to its upstream value. Sepa- 
ration on the external cowl surface occurs, extend- 
ing a significant distance upstream. 

The computations of Lai and Nelson^ were car- 
ried out for the nozzle experimental geometry of Re 
and Leavitt. The geometry is shown in Fig. 21, 
and the nozzle exhaust grids used are shown in 
Fig. 22. The spanwise grid distribution is shown 
in Fig. 23. A three-dimensional computation was 
carried out matching the nozzle operating condi- 
tions which exhaust into still air. Streamwise 
contours for Mach number are shown at four differ- 
ent spanwise positions in Fig. 24. Figure 24(a) 
shows the contours at the nozzle mid-plane, and 
Figs. 24(b) to (d) move progressively outward to 
the nozzle sidewall. Behavior of the shear layer 
formed on the three-dimensional surface of the jet 
is clearly evident, as one scans Figs. 24 and 25 
which are the cross-sectional planes in the x or 
downstream direction. The shear layer is observed 
to be highly three-dimensional in its structure 
with significant variations from top wall to bot- 
tom wall to side wall. Comparisons of the com- 
puted wall pressure with experimental data, 

Fig. 26, show good agreement for both surfaces. 

Improvements Required 

In all of the cases described so far, the 
nozzle entrance flow field was uniform and nonre- 
acting. In a scramjet module, the combustor and 
nozzle are very closely coupled which means that 
nonuniformities in temperature, pressure and spe- 
cies concentration will be present at the nozzle 
entrance. The sensitivity of the nozzle configura- 
tion to these nonuniformities has been a matter of 
debate for some time. It has been proposed by this 
author that three-dimensional reacting nonuniform 
profiles should be used in CFD codes to determine 
the performance sensitivity of the nozzle. In 
this way, it would be possible to determine what 
profiles are desirable at the nozzle entrance for 
high performance. The knowledge of those profiles 
would provide combustor designers with "targets" 
or goals to achieve at the combustor exit plane in 
terms of pressure, temperature and concentrations. 

A start on providing an answer to this approach 
has been initiated by Tsai and Yu^5 who performed 
three-dimensional Navi er-Stokes computations for a 
reacting hydrogen-air mixture as well as for frozen 
flow. Using a numerical scheme based on the LU 
approach of Yoon and Jameson, 36 pressure, tempera- 
ture and Mach number distributions for both flow 
solutions have been compared as well as the OH, 0 
and H mass fraction distributions. Calculation of 
the relative nozzle thrusts in a typical non axi- 
symmetric configuration is still required in order 
to support the proposal made by this author. The 
ability to compute multi-exhaust modules integra- 
ted into a realistic aftbody with reacting flow 
exhaust over a wide Mach number range has not yet 
been demons trated. 


Flow separation at transonic speeds will 
cause significant drag losses. To date, no realis- 
tic simulation of the transonic regime has been 
performed. Separation of the external cowl bound- 
ary layer and incorporation of means to prevent the 
separation also have not been simulated. Preven- 
tion methods could include burning or deployment 
of a wall section such as in the aftbody. As in 
the other components, turbulence and transition 
modeling remain as significant unknowns. Relami- 
narization of the aftbody boundary layer is also a 
phenomenon which could occur. Heat transfer to 
the aftbody is also a computational issue which 
requires a significant amount of attention. 

Conclusions 

Based on the results presented in this paper, 
it is evident that there are many portions or fea- 
tures of a hypersonic propulsion system that have 
not been properly or accurately simulated by numer- 
ical analyses. The inlet behavior and its perform- 
ance for both steady and unsteady operation have 
been computed for fairly simple geometries, and 
none of the computations incorporates all of the 
physics which are known to occur. The design of a 
high efficiency supersonic inlet, designed for a 
relatively narrow Mach number range, is not a sim- 
ple design matter. High recovery, reduced heat 
load, proper cowl lip design and low spillage drag 
are all required while delivering usable exit pro- 
file conditions to the combustor. In the spirit 
of a totally integrated propulsion system design, 
the features of the inlet flow field must be con- 
sidered as part of the combustor/nozzle design. 

Some thoughts on the utilization of inlet flow fea- 
tures for combustor design were presented in this 
paper. It is noted, finally, that unsteady inlet 
behavior with mass addition and flow spillage must 
be properly simulated by numerical methods in order 
to establish a reliable design. 

In the combustor region, it is clear that only 
crude design features have been modeled and com- 
puted. Realistic geometries, i ncorporati ng struts, 
wedges and wall injection must be computed in order 
to provide any design benefits. Aside from the 
additional geometric complexity needed, the funda- 
mental nature of turbulence-combustion interaction 
remains as a major unknown. The proper modeling 
of the unsteady flow field, both random and deter- 
ministic, with the combustion chemistry remains a 
major challenge. PDF modeling appears to be one 
of the promising avenues to pursue in this regard. 

Nozzle flows are characterized by intricate 
patterns of shock waves and shear layers. The 
asymmetry of the nozzle, combined with the shear 
layer position, reflect the shocks in various 
directions. Shear layer bending occurs at the 
shock intersection locations. The shape and struc- 
ture of these features vary considerably over the 
flight Mach number range with the attendant chang- 
es in nozzle back pressure. Much remaining work 
needs to be performed on the expansion of reacting 
gases with properly modeled entrance profiles of 
both flow stream aero properties and concentration 
profiles. It has been suggested in this paper that 
the efficient design of a hypersonic nozzle should 
proceed by exploring the best profiles for optimum 
nozzle performance. Those profiles then provide a 
goal for combustor exit conditions. Further defi- 
ciencies in the nozzle area are associated with 


7 



nonexistent modeling for reacting separated flow. 
Numerical orediction of the state of the boundary 
layer is also not feasible. In addition, accurate 
neat transfer computations for the nozzle wall 
have yet to be demonstrated over the range of oper- 
ating conditions. In closing, it is important to 
note that real gas and thermal and chemical non- 
equilibrium effects must be considered in all of 
the propulsion system components at the appropri- 
ate f 1 ight conditions . 

Advances in two other areas are required. The 
first advance involves the further development of 
accurate and efficient numerical methods required 
to improve solution accuracy with an attendant 
reduction in computing time. The second advance 
relates to computer technology and includes issues 
of speed, storage, structure and graphics. 

It is not the intent of this paper to present 
the opinion that a realistic computer simulation 
for a hypersonic vehicle poses an impossible task. 
Rather the intent is to assess, in as realistic a 
fashion as possible, what is achievable with todays 
knowledge, numerical codes and computers. The 
progress in these areas has been remarkable over 
the last decade and will continue to be so in the 
future. It is with this understanding that this 
author is confident that a complete simulation 
over the entire flight range will be possible some- 
day. However, the progression of the simulations 
to the point where design type information can be 
obtained, such as described earlier (Propulsion CFD 
Validation section), can only occur if the valida- 
tion process is carried forward with more realistic 
propulsion geometries. 
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Table 2. - Inlet Requi remencs 


Performance parameters 
Mass capture 

Kinetic energy efficiency 
Pressure recovery 
Heat load 

Exit profiles to combustor 

Computed variables required 
Wall pressure 
Skin friction 
Heat transfer 

Exit prof i 1e$ (u,v,w,T,Cj > 

Physical/chemical modeling requirements 
Shock boundary layer interactions 
Secondary/corner flows 
Mass injection 
Low Mach number bleed 
Flow separation 
Shock induced unsteadiness 
Turbulence, transition 
Equ i 1 1 br i urn chemi stry 
Flow unsteadiness 


Table 4. - Nozzle Requirements 


Performance parameters 
Thrust 
Moment 
Heat load 

Computed variables required 
Wall pressure 
Heat transfer 
Skin friction 

Physical/chemical modeling requirements 
Finite rate chemistry 
Turbu 1 ence 
Shock interactions 
Shear layers 
Secondary flows 
Separation 
Relami narization 
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FIGURE 14. - SIDEWALL PARTICLE TRACING, M = 12.25 (REF. i|). 



FIGURE 5. - ORIENTATION OF CROSSFLOW PLANE 150 LOCATED NEAR 
THE RAMP SHOULDER SHOWING PLANE OF SYMMETRY, SIDEWALL RAMP 
AND COWL SURFACES (REF. 12). 





1.0 CM 



FIGURE 10. - MACH NUMBER CONTOUR ON xy PLANE AT CENTER OF IN- 
JECTION PORT FOR CASE 1 (REF. 15). 



H 1 1 1 1 1 f 

3 CM 

FIGURE 11. - MACH NUMBER CONTOUR ON xy PLANE AT CENTER OF IN- 
JECTION PORT FOR CASE 2 (REF. 15). 
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FIGURE 13. - GRID SYSTEM PERSPECTIVE VIEW FOR 3-D COMPLETE GEOMETRY 
CALCULATION, CASE-1, a r = 20° (REF. 30). 
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FIGURE 14. - MACH CONTOUR PERSPECTIVE VIEW FOR 3-D COMPLETE GEOMETRY 
CALCULATION. CASE-1. a r = 20° (REF. 30). 



FIGURE 15. - SIMULATED SURFACE OIL FLOW FOR 3-D COMPLETE GEOMETRY 
CALCULATION. CASE-1, a r = 20° (REF. 30). 
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FIGURE 20. - DENSITY CONTOURS FOR CASE 3 (REF. 31). 



FIGURE 21. - NOZZLE GEOMETRY (REF. 32). 
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(d) I = 25, Z = 1.982. 

FIGURE 24.- STREAMWISE MACH NUMBER CONTOURS. CASE 1 
(REF. 32). 
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